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Abstract 

We find an exact solution for the stability limit of relativistic charged spheres 
for the case of constant gravitational mass density and constant charge den- 
sity. We argue that this provides an absolute stability limit for any relativistic 
charged sphere in which the gravitational mass density decreases with radius 
and the charge density increases with radius. We then provide a cruder absolute 
stability limit that applies to any charged sphere with a spherically symmetric 
mass and charge distribution. We give numerical results for all cases. In addi- 
tion, we discuss the example of a neutral sphere surrounded by a thin, charged 
shell. 
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1 Introduction 

Over the past decade, extremal black holes — black holes for which the charge equals 

the mass in geometric units — have been the subject of considerable interest, largely 

because such objects were the ones originally employed to derive the Bekenstein- 

Hawking entropy directly from string theory [1]. Aside from developments in string 
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theory, however, there has long been ample motivation to study extremal black holes 
because from the classical and semi-classical point of view they provide the "zero- 
temperature" limit in black hole thermodynamics. At the same time, substantial 
evidence suggests that one should not view extremal black holes as any sort of con- 
tinuous limit of their sub-extremal counterparts, black holes for which the charge is 
less than the mass. For example, the horizon structure of classical, charged black 
holes changes completely at extremality[2]. Some studies have also concluded that 
entropy is not well-defined for extremal black holes [3, 4]. More definitely, one knows 
from Israel's proof of the third law of black hole dynamics[5] that extremality cannot 
be attained in a finite time, and that the conclusion holds even under Hawking radia- 
tion and superradiance, which violate the assumptions of Israel's proof [3, 6] . Thus, it 
appears impossible to create an extremal black hole from a subcxtrcmal one, and the 
only remaining possibility is to produce one from the collapse of an already extremal 
object. 

For this reason it is of interest to investigate the stabihty of relativistic charged 
spheres. Previous studies along these lines have been mainly numerical [7, 6] and have 
concluded that while for Q < M collapse always takes place at a critical radius Rc 
outside the horizon, as Q approaches M, this critical radius approaches the horizon 
itself, The present paper is intended as an analytic companion to the numerical 
investigations. Our point of departure is the classic proof of Buchdahl[8], who showed 
that for uncharged stars gravitational collapse into a black hole will always take place 
when R < (9/4)M, regardless of equation of state. ■'^ (See Weinberg's Gravitation and 
Cosmology[10] for a clear presentation of Buchdahl's argument, or §3 of this paper.) 
Because in a charged sphere Coulomb repulsion tends to oppose the gravitational 
force, Rc should be less than (9/4) M. However, relativistically, the charge increases 
gravitational energy as well, and so at some point gravity always wins out and col- 
^ Throughout we use units in which G = c = 1. 
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lapse into a black hole takes place. At Q = M, the Coulomb repulsion equals the 
gravitational force and one finds numerically that Rc — R+. Thus one should have 
M <Rc< (9/4)M, always. 

Although for a given charge distribution one can indeed find Rc numerically, one 
does suspect that there must be an analytic proof, analogous to Buchdahl's, that ap- 
phes to relativistic charged spheres.^ In other words, given a value of Q/M we should 
be able to find an absolute bound on R/M, independent of other physical parame- 
ters, below which the object collapses into a black hole. Anninos and Rothman[6] 
(henceforth AR) intended to include such a proof as a supplement to their numerical 
investigation but as that project neared completion they learned that Yunqiang and 
Siming (henceforth YS) had already claimed to have given such a proof [9] . The YS 
demonstration, however, is far from transparent and does not appear to have ever 
been published. Moreover it does not provide a sharp value for the collapse radius, 
as in Buchdahl's proof, but rather gives a general lower bound on it. The important 
feature of this lower bound is that for Q < M it is always larger than as ex- 
pected from numerics. Not long ago we decided to take the opportunity to present 
a simplified version of this interesting result. In the process we have found an exact 
solution for the case of constant mass and charge densities, and this allows us, in a 
manner complete analogous to Buchdahl's, to put an absolute (sharp) stability limit 
on a very large class of objects, all those with charge density increasing radially, and 
gravitational mass density decreasing radially. This stability limit, which should cover 
essentially all cases of interest, is the main result of our paper. For the remaining 
cases, we present a proof similar to that of YS, but we hope with greater clarity, and 
give explicit numerical results for a lower bound on R^. We also give an exact solution 
for the stability limit of a neutral sphere surrounded by a charged shell. 

The paper is organized as follows. In Section 2 we introduce the relevant Einstein 

^In this paper we tend to speak of charged spheres rather than charged stars, as there is no good 
reason to think that charged stars, in the usual sense of the word, exist. 
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equations and introduce notation. In Section 3 we review the Q = case and derive 
an exact solution for the case of a neutral sphere surrounded by a charged shell. 
Section 4 is devoted to the main result of our paper: we solve exactly the case of 
constant charge density and constant gravitational mass density, derive its critical 
stability radius and show that it gives an absolute bound on the critical stability 
radii for all spherically symmetric distributions in which gravitational mass density 
gradient is negative and the charge density gradient is positive. In Section 5 we 
calculate a general lower bound on the critical stability radius. Finally in Section 6 
we summarize the results and draw conclusions. 

2 Einstein Equations 

As mentioned above, the plan is to find an absolute stability limit on R/M for rel- 
ativistic charged spheres that is independent of the equation of state and depends 
only on Q/M. We will assume throughout that the pressure p and density p are 
both positive, that the charge density is positive and that Q < M in order to avoid 
naked singularities; this last assumption ensures that spacetime is asymptotically 
predictable [2]. 

We also restrict attention to spherically symmetric mass and charge distributions, 
for which the metric can be written in the form 



where the metric components e^*^''^ and e^^*^*"^ are positive. 

As is well-known, the classic Reissner- Nordstrom (RN) solution for the charged 
spherically symmetric case gives 



where r > R and R is the outer radius of the sphere. The Schwarzschild-Droste (SD) 
solution is of course recovered by setting Q — 0. Both the RN and SD, however. 



(2.1) 




(2.2) 
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are vacuum solutions, concerned solely with the metric outside R. For the collapse 
problem we need to study the behavior of the metric functions A(r) and $(r) for 
r < R, where the pressure, the charge and mass densities are nonzero. The procedure 
for solving the "interior Reissner-Nordstrom equations" is nevertheless much the same 
as for the exterior case. One assumes (see AR or de Felice [7] for more details) a 
perfect-fluid stress-energy tensor for the hydrodynamic part 



hydro 



(2.3) 



while for the electromagnetic part 



MT^.)em = F;F,„ - \g,.F^pF'^^ = ^ diag[ e^*, -e^^, r^, r^sin^^ ]. (2.4) 

In these expressions p = Prm + e is the total mass density, prm is the rest mass density, 
e is the internal energy density, p is the fluid pressure, is the four- velocity and F^,^ 
is the electromagnetic field strength tensor. (T^u) em here is of the same form as for 
the exterior RN solution, as it must be by Gauss's law, except that instead of the 
total charge Q, we now have a q{r), the charge within any given radius r. Indeed, by 
definition 

q(r) = 47r T e^^''^+''^'V^f(r')dr', (2.5) 
Jo 

where is the charge density. (This is the usual definition of charge, modified only 
for metric curvature.) The boundary condition requires Q — q{R). 

With the above stress-energy tensor, the (00) Einstein equation is found to be 

g^(r) 



/2 



0<b' 
r 



p + 3p + 



47rr^ 



where / denotes derivatives with respect to r. Similarly, the (11) equation is 



p-p 



r 

and the (22) equation is 

.2A 1 + AV - $V = 47rr2e^^ p-p + 



q^{r) 



q^{r) 



(2.6) 



(2.7) 



(2.8) 
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The left-hand-side of these equations is necessarily the same as for the exterior SD 
or RN solutions; only the right-hand-side differs because of the nonzero stress-energy 
tensor. Following the standard procedure for deriving the exterior solutions, we can 
take linear combinations of Eqs. (2.6), (2.7) and (2.8) to eliminate the terms in 
One easily finds that for any r < R 

^-2A(.) ^ 1 _ 2m^ _ (2.9) 
r r 

where 

mAr) = in r p r'^ dr' and r(r) = f dr' . (2.10) 

Jo Jo r 

Here, mj(r) is the usual definition of the mass within a radius r. The subscript i 
denotes "internal" to emphasize that rUi contains both rest and internal energy. We 
use the designation because it will become necessary to distinguish mi{r) from the 
gravitational mass mg{r), defined momentarily. Requiring that (2.9) matches the 
exterior solution (2.2) at r = i?, gives 

or 

M = i/;(8.^^H-^)A-+g, (2.12) 

which defines the gravitational mass at it! (the mass measured by a satellite in orbit 
around the object). By Gauss's law, however, the same must be true at any radius, 
and so using the definition of T from Eq. (2.10), 

m,{r)=m,{r) + ^ + '^^, (2.13) 

which clarifies the distinction between and rUg. In terms of the gravitational mass, 
the metric function e^*^*"-* is 

=il- - =il- + . (2.14) 



Stability Limit ... 7 

One can write these functions either in terms of rrii and JF, or rUg and q, but because 
we do not in general know the charge distribution q{r) and hence J-{r), when thinking 
about boundary conditions it is much more convenient to use nig, since in that case 
e^^*"^ matches onto e^^^^ in the expected way. We will therefore generally use the 
second form. 

The pressure in these equations can be eliminated by taking three times Eq. (2.7) 
and adding it to Eq. (2.6), which yields 

r r \ J 
Upon multiplication by e~^^*/r, the left-hand-side turns out to be an exact differ- 
ential, and so, letting ({r) = e*^**^ as in Weinberg's notation [10], one has 



^2 rj^ ^5 



e\. (2.15) 



Eq. (2.15) will prove to be the fundamental equation of our analysis. It can be 
brought into perhaps more familiar form by noting by that 87rp/r = 2m\{r)/r^ and 
p[r) = q^r^. Then 



'mj(r)y^ 1 l'bq^{r) 3jF(r) 



C- (2.16) 



This is the equivalent of Weinberg's Eq.(11.6.14), employed in Buchdahl's proof of 
the absolute limit of stability for ordinary stars. Our equation, however, contains 
two more terms within the square brackets than the usual one, as well as an extra 
term in the expression (2.14) defining e^. We point out that although e^'^''^ has the 
simple form given in Eq. (2.14), no such closed form exists for C(r) for r < R. Indeed 
the differential equation (2.15) should be regarded as the equation defining ( in the 
interior of the sphere for given input distributions p and q. 

Requiring that ( and (' match on to the exterior RN solution at r = R gives the 
following important boundary conditions: 
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2M Q2 



) 



1/2 



( 



M g2 



) 



C{R) 



1 - 



R i?2 



(2.17) 



with M the gravitational mass in (2.12). 

Once C and $ have been computed using (2.15), one can easily show from (2.8), 
(2.9) and (2.13) that the pressure is given in terms of nig and q by 



Clearly, in order for the solution of (2.16) to be physically acceptable, we must require 
that C > and that p, as computed from (2.18), is nonnegative and satisfies a proper 
equation of state, usually assumed of the form p — p[prm]- For a given class of 
distributions p, q we define the critical stability radius Rc as the smallest possible 
radius for which a physically acceptable solution to (2.15) can be found. 

Then, in order to find a bound on the stability of the charged star, it is sufficient 
to show that if the radius is smaller than Rc, then no physically acceptable solution 
exists, whatever the choice of p, q. 



3 Neutral Sphere Surrounded by Charged Shell 



In this section we first review Buchdahl's Q = case, which illustrates the general 
strategy for finding Rc. The question then naturally arises as to whether the result 
changes for a neutral sphere surrounded by a charged shell. We show that the basic 
Q — scenario is fairly easily adapted to cover this case. 

For the moment, then, let us set Q = and drop the subscript i on m (there is 
now no distinction between between rUg and rrii). Eq. (2.16) then becomes 




(2.18) 




(3.1) 



where we have used e '^^^''^ = 1 — 2m{r)/r. 
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Let us now assume that for any physically reasonable star p' < 0. Then, because 
{m/r^y = 4Trr~^ Jq p'{x)x^dx < 0, it follows that 



in _ ^)i/2w 



Integrating this expression from r to R gives 



C')' < (3.2) 



1/ 2M^^/%,^^^ 1^^ 2mV/%, 



C'(i?)-i(l-^) C'(r)<0. (3.3) 



V R 

We now make use of the boundary conditions Eq. (2.17) with Q — 0. Inserting 
the expression for C'{R) into Eq. (3.3) gives 

C'M.|^(l-^)"". (3.4) 
Integrating again from to i? yields 

CiR) - C(0) > ^ / (3.5) 

SV J ^ - R3 Jo 2m(r)y/^ ^ ^ 

1 /2 

with ({R) = (^1 — . Now, in order to have ((0) > 0, we require that 



2My/2 M j-R rdr 



0<C(0)< 1-— --/ — (3.6) 



Note that, as remarked above, {m/r^)' < 0, and so m/r^ > M/R^ for all r < R. 
Plugging this into (3.6) we find: 



2M\^/2 1 

< 1 



(3.8) 



The integral is now trivially performed to get 

~R^) ~ 2 ~ V ~ ~R^) 
which immediately implies Buchdahl's result R> (9/4) M. 

Note that for stars with constant density all the above inequalities become equal- 
ities at the critical radius, and so the value (9/4)M is precisely their critical stability 
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radius.3 The value = (9/4)M gives an absolute (sharp!) stability limit for all stars 
with distributions satisfying p' < 0. If any such star is compressed to the point that 
R < Rc gravitational collapse necessarily takes place. 



Let us now modify the above computation to handle the case of a neutral sphere 
of constant density surrounded by a thin shell of internal ( "inertial" ) mass K that 
carries a uniformly distributed charge Q < M. In such a situation, Eq. (2.14) shows 
that A suffers a discontinuity at r = R. Let R~ and i?"*" represent the inner and outer 
radii of the shell, Mj„j be the mass interior to the shell and Mg = K + Q'^/2R be 
the gravitational mass of the shell (cf. Eq. (2.12)). Then M = Mj„t + Mg is the 
total mass and e-^^^-^") = 1 - 2Mint/R, while e-^^^^^) = 1 - 2M/R + Q'^/R'^. Since 
the jump in A is finite, however, Eq. (2.18) implies that any discontinuity in $' and 
hence in (' is finite as well. Thus ( itself is continuous at the boundary with precisely 
the value given by the first of Eqs. (2.17). 

One can greatly simplify the calculations by assuming that K — 0, in which case 
A and $ are both continuous oXr = R and only their derivatives suffer a discontinuity 
at the surface. (For a more detailed discussion of these issues we refer the reader to 
Cohen and Cohen [11], who derive the solution for a thin charged shell of radius R, 
with Mint = 0.) 

So, let us take K = and assume that interior to the shell p = 3 Mint / {4:71 R'^), 

with Mint — M — Q'^/{2R). In order for the shell to be stable against gravitational 

collapse, it is necessary to have a nonzero elastic stress tensor concentrated on the 

surface, as assumed in [11]. This means that the stress energy tensor must be modified 

by the addition of a term (T^iy)ez whose only nonzero elements are Tg0 = r'^S6{r — R) 

and T^<^ = sin^ 9SS{r — R), where S is the elastic energy and the delta function 

■^It is straightforward to check that in the uncharged case, if C > 0, then p > and finite, Vr < R. 
This means that the only condition to be imposed for the solution to be physically acceptable is 
C > 0, in other words, precisely the condition we imposed above. 
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is normalized such that / dr47rr'^6{r — R) = 1. The presence of {T^y)ei modifies the 
Einstein equations as follows: Eq. (2.6) contains an extra term 87re^^5'(5(r — R) on 
the right hand side; Eq. (2.7) contains an extra term — 87re^^5'5(r — R) on the right 
hand side; Eq. (2.8) is unchanged. With these additions Eq. (2.16) becomes 

mj(r) Y ^ 1 ( hq^{r 



r / 



A(r) 



1 

^ 2^ 



r r 



where 



-2A(r-) 



1 — 2mi{r)/r r < R 

1 - 2M/r + gVr^ r>R 



C, (3.9) 



(3.10) 



and C = e ^ for r > i?. Integrating both sides of Eq. (3.9) between R and 
R^, and using the fact that A and $ are both continuous at r = i?, we see that 
e-^(^)(C'(i?+) - C'(^")) = Ke^S/R\ or 



^1 - 2Mint/R L 

Computing C,'{R~) from Eq. (2.18) yields 

yjl - 2MijR 

and comparing the two expressions shows that 



int 



i?2 2i?3 



25 



(3.11) 



(3.12) 



inRp- 



2S 



(3.13) 



2i?3 i?2 

The parameter 5" should be chosen such that p_ > 0. 

We can now proceed as in the Q — case, integrating (^^e~^('^ — from 
to R~ with the new boundary conditions ({R 
+ ^nRp 



l-2Mint/R and ('{R') = 
The result is that the new critical radius is smallest when 



p_ = 0, in which case 5" = —Q^/AR, consistent with the result in [11] for K — {). We 
then find that Rc is precisely (9/4)Mj„t, as one might expect from Gauss's law. With 
Mint = M — I (2i?c); solving for R^ in terms of M and Q gives 



(3.14) 
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For Q = 0, i?c = 9/4M, as expected, and a nonzero Q indeed lowers Rc- The extremal 
case, Q — M, gives Rc — 3/2M, which is plausible, as below we will find that for the 
full extremal charged sphere Rc — M. 

4 Constant Density Case 

The special case of perhaps greatest interest (and the easiest one to handle), is that 
of constant density, nig oc and q oa r^. Remarkably, we are able to find an 
exact solution in this situation. Moreover, because a neutral test particle senses the 
gravitational mass nig within a radius r, it is evidently nig that plays the role did 
in the Q — case. In other words, it is nig that determines the weight of material in 
the sphere, and a physically reasonably requirement for stability is that < 0, where 
pg is the gravitational mass density. If we additionally impose the requirement that 
the charge density is positive and increases outwards (that is, q' >^ and (q/r^)' > 0), 
which also seems reasonable if like charges repel, then we will also be able to find, in 
complete analogy with Q = case, an absolute lower bound on the critical radius of 
any charged sphere meeting the two conditions. 

We begin by rewriting the fundamental equation (2.16) in terms of m^: 



2{^.\\ -1/2 



nig_\ _jq_ 

^3 / V ^4 



C (4.1) 



With the ansatz that nig — M{r/R)^ and q — Q{r/RY, the first term in the square 
brackets vanishes and Eq. (4.1) becomes 



where now 

P-^- = w 1 - 9 , , , 

R \RJ i?2 \R 

Here and in what follows the subscript c refers to "constant-density case." 



e-=>l-2^f^)M;f^V. (4.3) 
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Let us define a new variable ( such that 

C(/c(r)) = C(r) ; Ur) ^ R-^ f dx xe^-^^^ 

Jo 

Substituting (4.4) into Eq. (4.2) gives at once 
which has the obvious solution 



Cc(/c)=Cie'5/=/'' + C2e 



-Qfc/R 



Moreover, /c(r) is a standard integral: 

x/R^ 



13 



(4.4) 



(4.5) 



(4.6) 



fc{r) = / dx 



l-2f (*) +f (I) 



4 2M 



dy 



R I M Q I ^ 



Vl-2y + (^)y2 



-2g^°MQ--^y + yi-2y + (^)y^ 



(4.7) 



or 



/c(^) = 7^ log 



M/Q + 1 



2g M/g - gr2 /i?^ + e-^r) ■ 



(4.8) 



Imposing the boundary conditions Cc{fc{R)) — e ^''^^^ and d(^c{fc{R))/dfc 
M/R — Q'^/R'^ we find after some algebra 



i(M/g-g/i? + e-^=(«))3/2 



C2 



M/g + 1 
= -^(M/g-g/i?-e-^=(^) 



M/g + i 



\ M/Q-Q/R + e-^'^i^y 



(4.9) 



and so, finally, the exact solution for C is 

1 r(M/g-g/i? + e-^^(«)) 



3/2 



Cc(r) 



'^[^M/Q-QryR^ + e-^^ 



(r) 
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The condition for this solution to be physical is Cc(0) > 0.^ As in the Q — case we 
get the equation for the critical radius by setting Cc(0) = 0, which yields: 



One easily sees from Eq. (4.11) that Q = M imphes Rc = M, as claimed in the 
Introduction. For other values of Q/M we solve this equation for R/M. Letting 
II — M/R and a — Q/M in Eq. (4.11), we find after some further algebra: 



Thus the exact solution for the critical radius boils down to finding the roots of this 
cubic equation for /i in terms of a. One immediately sees that Q — implies that 
Rc/M — 9/4. Numerical results for various values oi Q/M are given in Table 1 and 
plotted in Figure 1. 



straightforward computation shows that if Cc > then p, as computed from (2.18), is auto- 
matically positive and finite, as expected. 



(M/g - Q/R + e 



-Ac(R) 



)' = (M/g + 1) (M/g - Q/R - e-^=(^)) (4.11) 



4(7V^ - 12(jV^ + (9 + 3(j2)/x -4 = 0. 



(4.12) 
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0/M 


Rc/M 





2.250 


.1 


2.244 


.2 


2.226 


.3 


2.196 


.4 


2.152 


.5 


2.093 


.6 


2.016 


.7 


1.915 


.8 


1.781 


.9 


1.586 


.99 


1.224 


.999 


1.091 


.9999 


1.039 



Table 1. The stability limit Rc/M for the const ant- density sphere, tabulated for 

various values oi Q/M. 




0.2 0.4 0.6 0.8 1 
Q/M 



Figure 1: The stability limit Rc/M plotted against Q/M for the constant density case. 

Thus, we have found an exact solution for the constant density case. We now ar- 
gue, as indicated earlier, that any physically reasonable solution should have p'g < 
and, for a positive charge density, {q/r^)' > 0. If so, the constant density case 
maximizes the expression in square brackets in the right-hand-side of Eq. (4.1) 
and then ((l/r)e-^C')' < e^{Q'^r/R^)C, as long as C > 0. In terms of C and 
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f — R'"^ Jq dx xe^^^\ Eq. (4.5) is now replaced by the inequality 

< ^C. (4.13) 
dp - R'^ ^ ' 

Note that now / ^ fc, because the e'^ appearing in the definition of / is no longer that 
of the constant-density solution (4.3) but instead the general (1 — 2mg/r + /r"^) ^^"^ . 

As discussed in Appendix 1, the differential inequahty (4.13) implies that C < Co, 
for all / < /(-R), where Co(/) is the solution to the differential equation d'^Co/dp — 
{Q'^/R'^)(o satisfying the same boundary conditions as ( does. One finds that Co(/) = 
cOgQ//^ + cOe-Q^/«, with 



<-2 — 



2 \Q R 

Note that the conditions p'g ~ (nig/r'^y < and (q/r^)' > imply that e^^^*^'-' is 
always smaller than e-2^= = {l-2Mr^ / R^+Q'^r^/ R^). Therefore Co(/(r)) < Cc(/c(r)). 

The conclusion is that C{f{r)) < Co{f{r)) < Cc(/c(r)) or, equivalently ({r) < Qir), 
where ({r) is the general solution to (4.1). Consequently, the critical radius for any 
distribution with < 0, g' > and (q/r^)' > is always larger then the critical 
stability radius plotted in Fig.l. The constant density case therefore provides us 
with an absolute stability limit for any relativistic charged sphere satisfying these 
conditions. 

Of course, we do not know exactly what a "realistic" charge distribution is for 
such objects, and so below we provide a somewhat cruder bound that is independent 
of any assumptions whatsoever about the charge distribution. The remainder of the 
paper consists of a formal proof along the lines of YS, but we have streamlined the 
presentation, filled in a number of gaps and present explicit numerical results. 
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The plan is now to bound the behavior of the solution to Eq.(4.1) for < Q < 
M under the most general conditions possible. To reiterate, we assume that any 
physically acceptable solution meets only the following conditions: 

p(r) > , C{r) > , 

Q<Q <M , R> R+ . (5.1) 

Here i?+ = M+i/M^ — Q'^ is the outer horizon of a RN black hole; if i? = i?+ then 
gravitational collapse of the charged sphere has already taken place. The quantities 
rrii, nig and q will be considered inputs that are related by Eq.(2.13) and they satisfy 
the conditions 

mg{R) = M , q{R) = Q , 

Trig > q , Trig + rrig — q^ < r . (5.2) 

The last two conditions are required to avoid naked singularities, as discussed in 
Section 2. 

Once again, we define the critical instability radius Rc{M, Q) as the smallest radius 
R > R+ for which a physically acceptable solution can be found in [0, R] for any input 
functions mg,q satisfying (5.2). We also parameterize the difference between Q and 

M by 

A = ^1 - g7M2 < 1 (5.3) 

and assume R = R+{1 + e), where e is some number (not necessarily small). With 
this notation, R^ — (1 + A)M and 

R= {l + {l + e)A)M. (5.4) 

We shall prove below that, under the conditions (5.1)-(5.2), the critical stability radius 
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admits a general lower bound of the form 
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Rc>{l + {l + So)A)M (5.5) 

with Eo ^ 1/264. We do not believe that this specific value of Sq has any physical 
relevance: it is merely a byproduct of our estimates and it can certainly be improved. 

Let us now turn to the proof of (5.5). The strategy will consist in demonstrating 
that, whenever i? < (1 + (1 + eo)A.) M, for some suitable Eq to be constructed below, 
then every possible solution ( of Eq. (4.1) becomes negative somewhere in [0, R), for 
any rrig, q satisfying conditions (5.2). Thus the solution becomes physically unaccept- 
able. In the proof we shall need the following preliminary estimate: 



Lemma 1 //i? < (1 + (1 + So)A) M and 6= {el + 2£o)^/^ then 



For a proof of Lemma 1, see Appendix 2. 

We now begin to study Eq.(4.1). Integrating both sides between r and R and 
using the boundary conditions (2.17) we find, in analogy to Eq. (3.3): 



R^\R i?2 



-A(r) 



' — q\x) ' 



x^ 



Cix) (5.7) 



A crucial point in our proof consists in finding a uniform bound on the right-hand- 
side, independent of (, rUg and g, at least for r close enough to R. This will allow us 
to dispense with the details oi C,. In Appendix 3 we prove the following key estimate. 



Lemma 2 Let (5 < I. If R<{1 + {1 + £o)A) M and Sq = /3^q;o, with 

(^-1)V4 



1 + ^1 + (V2-1)V4' 
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then either the r.h.s. of (5.7) is uniformly bounded from above by — for 

all r e [/3R, R] or C(r) = for some r e [/3R, R). 



Remark. Note that if Sq is chosen as in Lemma 2 then the constant 5 in Lemma 
1 satisfies 6 < /5^(\/2 — 1)/2 (see the proof of Lemma 2). This fact will be used below. 

Now assume that C > for any r < R. By Lemma 2, for any r e [PR, R] with 
(3 < 1, if i? < (! + (! + £o)A) M where now Eq = /3^ao, then the right hand side of 
(5.7) is uniformly bounded above by — As a consequence, from (5.7) we 

find 

C'(r) >!:^(M_Q^\ 

' - 2R'^ \R Ry' 
Integrating this inequality between f^R and R and using > 1, we get: 



+ ^-C(/?i?)>2^-^W^, (5.8) 



R i?2 ^^"^ ' - 2\R R'^ J JpR i?2 ' 
or 

0<C(/Jfl)<^l-- + y;-^(^--^j (5.9) 

On the other hand, under the assumption that R < (1 + (1 + 80)^) M with Eq ~ 
P'^ao, we have from Lemma 1 and the remark after Lemma 2, 



If we merely set (3'^^^^^ — which gives 



we then have a contradiction, and this implies that no physical solution to the Einstein 
equations can be found. This proves (5.5), with an explicit bound on Eq given by /5qQ;o; 
or £0 — 1/264. 



Stability Limit ... 

We can now tabulate Rc/M 
results are shown in Table 2. 



20 

> 1 + (1 + £o)^ for various values of Q/M. The 



Q/M 


Rc/M 





2.250 


.1 


1.999 


.2 


1.983 


.3 


1.956 


.4 


1.920 


.5 


1.869 


.6 


1.830 


.7 


1.717 


.8 


1.602 


.9 


1.437 


.99 


1.142 


.999 


1.045 


.9999 


1.014 



Table 2. The lower bound on Rc/M is tabulated for and various values of Q/M. 

Because these figures represent a lower bound on Rc/M, they should all lie beneath 
the corresponding numbers of Tabic 1, and indeed they do. Any relativistic charged 
sphere, regardless of equation of state, must have a critical stability value of Rc/M 
greater than the values presented here. The discussion above fails for Q = M, in 
which case the stability bound is simply R > R+. 

6 Conclusions 

The main result of this paper is an exact solution for the stability limit of constant- 
density relativistic charged spheres for all Q < M. We also argued that in any 
"physically reasonable" case where the gravitational mass density decreases with the 
radius and the charge density increases with radius, the constant-density case provides 
an absolute stability limit for all charged spheres. If calculation of a stability limit for 
rotating objects proves tractable we would expect a quantitatively similar behavior 
for Rc/M, given a value of a/M, the angular-momentum-parameter-to-mass ratio. 
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For the most general charged case we found a cruder bound that is independent 
of any assumption about the mass and charge distribution, except for the basic con- 
ditions (5.2). Both "physically reasonable" and general bounds approach the horizon 
i?+ in the limit Q M. That in this limit the critical "collapse" radius is precisely 
i?+ is intuitively reasonable because Q = M is the point at which the Coulomb repul- 
sion equals the gravitational force. Such a state is evidently unstable, but our results 
apparently do not exclude what seems to be the only route to producing an extremal 
black hole: to first create an extremal charged sphere, then compress it to the hori- 
zon. Nevertheless, given the other issues surrounding extremal bodies mentioned in 
the Introduction, one should continue to hesitate before regarding such objects as a 
smooth limit of the sub-extremal state. 

Note Added: Just prior to submission of this paper we learned that a similar proof 
of our main result has been independently given in [12]. 

Acknowledgments. The work of A.G. was partially supported by U.S. National 
Science Foundation grant PHY 01 39984, which is gratefully acknowledged. 

Appendix 1: On the differential inequality (4.13) 

In this Appendix we prove that if satisfies inequality (4.13) with boundary 
conditions C(/(i?)) = e-^(«) and dCifm/df = M/R - / R\ then ((/) < Co(/), 
for all < / < /(-R). Again, Co(/) is the solution to the differential equation 
(PCo{f)/dp — {Q^ / R^)Co{f) satisfying the same boundary conditions as C(/) does. It 
follows that the constant density case nig ex r^,q oc provides an absolute stability 
limit for any relativistic charged sphere in which < 0, g' > and (q/r^)' > 0. 
In order to prove that ( < (q, we define g{x) = Co(/(-R) — x) — ({f{R) — x) and 
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show that g{x) > 0, for all < a; < f{R). Note that g{0) = dg{0)/dx = and from 
(4.13) that cPg/dx^ > {Q'^/B?) g, for all < x < f{R). Without loss of generahty, 
we can assume that g{x) is not identically zero in a right- neighborhood of the origin. 
Let us pick some xq < R/Q (the reason for this choice will become clearer below) and 



(1) g{xi) — 0. This means that g{x) must be nonnegative in [0, xq], which in turn 
implies that g{x) is nonnegative for all < x < f{R). If this were not the case, 
then there would be some x*, for < x* < f{R), such that g{x*) > would be a 
local maximum. But this would mean that (P g{x*) / dx"^ < 0, which contradicts the 
condition d'^g{x*)/dx'^ > (Q^/R^) g{x*) > 0. 

(2) g{xi) < 0. By the mean-value theorem g{xi) = Xidg{x2)/dx, for some < 
X2 < Xi. Again by the mean- value theorem, dg{x2)/dx — X2(P' g{x^) / dx"^ , for some 
< X3 < X2. On the other hand d'^ g{xz) / dx'^ > {Q'^ / R^) g{x^), and so we have 
> g{xi) > XiX2{Q^ / R^) gixs) , implying in particular that g{x3) < 0. Since g{xi) is 
the minimum of g{x) in [0, Xq] we also have |5'(a;3)| < |g'(xi)|. Finally: 



Because < X2 < Xi < Xq < R/Q, we see that the r.h.s. of this inequality is strictly 
smaller than |g'(xi)|, but this is a contradiction. Thus g{xi) < is an impossibility 
and the proof is concluded. 



Appendix 2: Proof of Lemma 1 

Let R^{1 + {1 + e)A)M as in Eq. (5.4) for < £ < Sq. Then in terms of A and s: 



let < Xi < Xq satisfy g{xi) = min^-gjo^a;!,] g{x). There are two cases: 



0<\g{xi)\<xiX2{QyR^M^i)\ 



(Al.l) 



M 1 



R ! + (! + £) A 



g2 _ 1- A=^ 

R^~ [! + (! + £) A]' 



(A2.1) 
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Inequality (5.6), which we want to prove, now takes the form: 



1 ^+ .,<S^( ' X. (A2.2) 

1 + (1 + £)A [l + {l + e)Af- \l + {l + e)A [1 + (1 + £)A]V 

Multiplying both sides by [1 + (1 + s) A]^ gives 

Notice that for < Q < M the right-hand-side is always > S^. So if we choose 
6'^ — el + 2£o, Eq. (A2.3) is certainly satisfied and the lemma is proven. ■ 

Appendix 3: Proof of Lemma 2 

Let us denote the right-hand-side of (5.7) by G{r). The integral vanishes if the lower 
limit is R; hence G{R) = 0. So, by continuity, G{r) will be less than ^ — in 
a small enough interval of the form [ro, R]. Let us pick some (3 <1 and let Sq — P^ao, 

with an = /'^ ^ (^^jig relevance of this specific number will be made clear 

below). We want to show that either G{r) < 2W ~ ^) ^ 
({r) = for some r e {/3R,R), in which case the solution becomes unphysical. 

We proceed by contradiction. Assume that ({r) > in [PR,R] and that G{r) — 
2^ ~ w) changes sign in the same interval. This means that there is some P' > (3 
such that both GiP'R) = 2I2 (f - i^) and G{r) < 2^ (f - 1^), Vr e {f3'R,R]. 
Now, from the definition of (Eq. (2.9)), one has 



dx 



(A3.1) 



so that G{j3'R) can be immediately rewritten as: 



G(f3'R) = r dx 

J,3'R 



1 d{e 



Cix) + dx e^(^) f-2^ + C(.x-) (A3.2) 



I3'R \ X^ X^ 



x"^ dx 

By (5.2), using in particular q{x) < mg{x) < x, we find that —2mg/x^ -\- 3q^/x^ < 
1/x^. Thus, after integrating the first term by parts. 



Gif3'R) < -^^e-^(^) + ^e-^(^'^) + dx e'^ f ^ - 2^] + ^ dx'- 

R^ [p Ry Jp'R \X-^ X-^ J J/3'R 



X^ 



(A3.3) 
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Note that from (2.18), under the condition that p > and r > rrig > q, we must have 
C > for C > 0. So, neglecting the negative terms in (A3. 3), and using the fact that 
< 1, C' ^ ^iid that x > P'R, the inequahty becomes: 

Now, with the supposition that G((3'R) — — as well as the boundary 

condition (2.17) for ^, 

dx e^(^) > ^^~f^Q^ - 2/5'^ (A3-5) 

On the other hand, in view of (5.7) and of the condition that G{r) < — 
Vr e {P'R,R], we have 

Integrating between (3'R and i? yields 

and employing the boundary conditions on C(-R) once again gives 



Since inequality (A3. 5) gives a minimum for the integral in this expression we can 
insert (A3. 5) into into the right-hand-side here to get, finally. 



1 2M , 

^ R ^ B? 



2(3'R . 
(A3.9) 



If we let y — _ 2M _|_ _ then the previous expression becomes 



y' + {l3'ry-l{P'r>0, (A3.10) 
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which imphes that 

y>ir^>f^. (A3.n) 

On the other hand, we see that Lemma 1 states that y < S, where 5 = {el + 2eo)^''^- 
Take Eq smaller than the positive root of the equation x"^ + 2x — (3^{\/2 — 1)^/4 = 0, 
that is, smaller than x+ = -1 + ^^1 + /^^(^ - 1)2/4. For example, Eq = /3^/47.13 < 
- 1) V[4 + 2^4 + (V2- 1)2] does the job. Then y < (3'^^. We have reached 
a contradiction, and the Lemma is proved. ■ 
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